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Systematic  series  of  hulls  based  on  a  parent  hull  have  long  been  used  by  naval  architects  and 
hydrodynamicists  to  determine  the  effects  of  hull  form  on  flow  and  resistance.  This  memoran¬ 
dum  describes  an  interactive  computer  program,  MODHLL,  which  creates  modified  hulls  from 
an  existing  parent  hull.  Although  the  parent  hull  requires  a  large  data  file  for  its  definition, 
each  modified  hull  is  defined  by  only  a  few  extra  words  of  data  since  only  the  way  it  differs 
from  the  parent  need  be  defined.  Thus  a  whole  series  of  hull  forms  may  be  defined  using  only 
a  little  extra  storage  space  than  the  parent. 

The  program  is  based  on  a  generalization  of  Lackenby’s  method  in  which  the  hull  is 
modified  by  moving  the  stations  between  the  fore  and  aft  perpendiculars.  Each  modified  hull 
is  specified  by  choosing  values  for  the  length,  breadth,  draft,  depth  and  either 

'  1.  the  fore  and  the  aft  prismatic  coefficients,  or 
2,  the  total  prismatic  coefficient  and  the  longitudinal  centre  of  buoyancy. 
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Un  ensemble  systematique  de  coques  basees  sur  une  coque  mere  est 
dupuis  longtemps  utilise  par  les  ingenieurs  des  constructions  navales  et  les 
specialistes  en  hydrodynamique  pour  determiner  les  effets  de  la  forme  des 
coques  sur  i'ecoulement  et  la  resistance.  Ce  memoire  decrit  un  programme 
d'ordinateur  interactif,  le  MODHLL,  qui  produit  des  formes  de  coques  modifiees 
a  partir  d’une  coque  mere  exist’ante.  Quoique  un  imposant  fichier  de  donnees 
soit  necessaire  pour  la  definition  de  la  coque  mere,  chacune  des  coques 
modifiees  n'est  definie  qu'au  moyen  de  quelques  mots  supplementaires  puisqu'il 
n’est  necessaire  de  definir  que  la  maniere  dont  elles  different  de  la  coque 
mere.  Ainsi  il  est  possible  de  definir  tout  un  ensemble  de  formes  de  coques 
en  n'utilisant  qu'un  peu  plus  d'espace  de  stockage  qu'il  n'en  faut  pour  la 
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coque  mere. 

Le  programme  est  base  sur  une  generalisation  de  la  methode  de 
Lackenby  par  laquelle  la  coque  est  modifiee  en  depla?ant  les  couples  de  trace 
entre  les  perpendiculaires  avant  et  arriere.  Chacune  des  coques  modifiees  est 
specif iee  par  le  choix  de  valeurs  pour  la  longueur,  la  largeur,  le  tirant 
d’eau,  la  hauteur  et 

1.  les  coefficients  prismatiques  avant  et  arriere  ou 

2.  le  coefficient  prismatique  total  et  le  centre  de  carene  longitudinal. 
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Notation 


a,b,c,...,k,l : 
A(x): 

B: 

Bi(z): 

c,-: 

D: 

<U: 

e.: 

/*(*>  y,...): 
9(x,y,..-): 
g*[x,y,...): 
k: 
L: 
N: 
P\,Pl- 
Pi,  PI¬ 
S'. 

sgn(i): 
U  • 


coefficients  of  the  polynomial  segments  defining  the  hull  transformation 

the  sectional  area  at  station  x 

the  half-breadth  of  the  hull 

the  ith  B-spline  basis  function  of  order  4 

the  ith  B-spline  basis  function  of  order  k 

weights  used  in  Gaussian  quadrature 

the  depth  of  the  hull 

the  derivative  of  the  function  at  poin  x,- 

the  error  associated  with  the  ith  data  point  of  a  spline  function 
a  function  of  x,  y, . . . 

first  partial  derivative  of  /(x,  y, ...)  with  respect  to  x:  d/(x,  y, . .  .)/dx 
a  function  of  x,  y, . . . 

the  first  partial  derivative  of  y(x,  y, . . .)  with  respect  to  x:  dg(x,  y, . .  .)/dx 
the  order  of  a  spline 

the  length  of  the  hull  between  perpendiculars 

the  number  of  B-splines  used  to  generate  a  spline  curve 

the  breakpoints  of  the  transformation  x'(x) 

the  stations  on  the  original  hull  which  correspond  to  p\  and  pi  respectively 
non-cartesian  hull  coordinate 
sign  of  t 

the  roots  of  the  Legendre  polynomial  of  fourth  order 
function  defining  the  waterline 
station  number 


X:  cartesian  hull  coordinate 
Xg:  longitudinal  centre  of  buoyancy 
Y :  cartesian  hull  coordinate 

y(x,  s):  function  defining  the  Y  coordinate  of  the  hull  in  terms  of  the  non-cartesian 
coordinates  (x,  s). 

Z:  cartesian  hull  coordinate 

2(1,5):  function  defining  the  Z  coordinate  of  the  hull  in  terms  of  the  non-cartesian 
coordinates  (x,  3). 

Z4:  the  design  draft  of  the  hull 

Zw:  the  draft  of  the  hull 

a,-,  7,-,  A,-:  parameters  used  in  the  test  for  monotonicity  of  the  transformation  x'(x) 
ft:  spline  coefficients  for  the  waterline  curve  or  curve  of  sectional  areas 
aft  prismatic  coefficient 
<j>F‘.  forward  prismatic  coefficient 
<f>:  total  prismatic  coefficient 
X 2 :  the  chi-square  of  a  spline 
spline  coefficients  (i  =1,4) 

Note:  Variables  with  primes  refer  to  the  parent  hull:  e.g.  X'B  is  the  longitudinal  centre  of 
buoyancy  of  the  original  hull,  while  Xg  is  the  longitudinal  centre  of  buoyancy  for  the  modified 
hull. 


1  Introduction 


Systematic  series  of  hulls  based  on  a  parent  hull  have  long  been  used  by  naval  architects 
and  hydrodynamicists  to  determine  the  effects  of  hull  form  on  flow  and  resistance.  Since 
computer  programs  are  now  used  for  calculating  detailed  hydrodynamic  performance  of  hulls, 
an  efficient  means  of  generating  computer  representations  of  series  of  hulls  is  both  useful  and 
convenient  for  determining  the  relations  between  hull  form  and  performance.  The  suite  of 
programs  HLLFLO[l],  developed  at  DREA,  may  be  used  to  calculate  the  flow  around  a  ship 
hull;  however,  HLLFLO’s  original  system  for  representing  hulls  in  the  computer[2,3]  required 
the  generation  of  a  separate  representation  of  each  hull.  This  is  not  only  a  time  consuming  task, 
but  also  can  require  large  amounts  of  disk  storage  since  each  hull  is  stored  in  a  file  containing 
several  thousand  words  of  data. 

This  memorandum  describes  an  interactive  computer  program,  MODHLL,  which  creates 
modified  hulls  from  an  existing  parent  hull.  Each  modified  hull  is  defined  by  only  35  to  70 
words  of  data  (to  be  compared  with  the  several  thousand  for  the  parent)  since  only  the  way  it 
differs  from  the  parent  need  be  defined.  Thus  a  whole  series  of  hull  forms  may  be  defined  using 
only  a  little  extra  storage  space  than  the  parent. 

MODHLL  is  based  on  a  generalization  of  Lackenby’s  method[4]  in  which  a  hull  is  modified 
by  moving  the  stations  between  the  fore  and  aft  perpendiculars.  Unlike  Lackenby’s  method, 
MODHLL  also  allows  the  draft  to  be  changed  independently  of  the  hull  depth.  Each  modified 
hull  is  specified  by  choosing  values  for  the  length,  breadth,  draft,  depth  and  either 

1.  the  fore  and  the  aft  prismatic  coefficients,  or 

2.  the  total  prismatic  coefficient  and  the  longitudinal  centre  of  buoyancy. 


Modification  of  the  draft  leads  to  some  confusion  of  terminology  which  is  best  dealt  with 
at  the  outset.  On  the  parent  hull,  the  “draft”  and  the  “design  draft”  are,  by  definition,  one  and 
the  same;  however,  on  a  modified  hull  the  “draft”  and  the  “design  draft”  may  be  different.  The 
draft  of  a  modified  hull  is  a  value  specified  by  the  MODHLL  user  and  may  be  used  to  account 
for  changes  in  hydrostatic  condition.  On  the  other  hand,  the  design  draft  of  a  modified  hull  is 
defined  to  be  the  design  draft  of  the  parent  scaled  by  the  relative  change  in  the  depth  of  the 
two  hulls.  That  is, 


Zi 


Z'dD 

D' 


(1.1) 


where  Zd  is  the  design  draft  of  the  modified  hull,  Z'd  is  the  design  draft  of  the  parent,  D  is  the 
depth  of  the  modified  hull,  and  D1  is  the  depth  of  the  parent.  The  “design  waterline”  is  the 
waterline  which  yields  the  design  draft;  the  “draft  waterline”  is  the  waterline  which  yields  the 
draft. 
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All  HLLFLO  programs  will  accept  the  modified  hull  representations  generated  by 
MODHLL  in  lieu  of  the  parent  file;  whether  a  hull  is  modified  or  is  a  parent  is  irrelevant 
to  the  remainder  of  the  HLLFLO  programs. 


2  Coordinate  Systems 


The  hull  is  embedded  in  a  cartesian  coordinate  system  whose  axes  are  aligned  with  the 
waterlines  and  centreplane.  The  three  cartesian  coordinates  are  denoted  ( X ,  Y,  Z)  (see  Figure  1) 
and  are  defined  as  follows. 

•  X  increases  from  bow  to  stern.  X  =  0  is  the  forward  perpendicular  and  X  =  L  is  the  aft 
perpendicular,  where  L  is  the  length  between  perpendiculars.  The  planes  of  constant  X 
are  hull  sections. 

•  Y  increases  from  the  centreplane  to  the  point  of  maximum  breadth.  Y  =  0  is  the  centre- 
plane.  Owing  to  the  symmetry  of  the  hull,  it  is  usually  not  necessary  to  decide  whether  Y 
increases  to  port  or  to  starboard;  however,  to  maintain  a  right-handed  coordinate  system, 
Y  should  increase  to  starboard.  Surfaces  of  constant  Y  are  parallel  to  the  centreplane. 

•  Z  increases  from  the  baseline  to  the  point  of  maximum  height  above  the  baseline.  The 
surfaces  of  constant  Z  are  waterplanes.  Z  =  0  is  the  waterplane  which  contains  the 
baseline. 

The  computer  representation  of  the  hull  requires  a  two  dimensional  coordinate  system 
(x,  a)  which  specifies  points  lying  on  the  surface  of  the  hull.  Given  a  point  (x,s)  on  the  hull 
surface,  one  must  be  able  to  calculate  the  cartesian  components  (A,  Y,  Z)  of  that  point.  X,Y, 
and  Z  can  then  be  regarded  simply  as  functions  of  x  and  a.  In  the  system  currently  used  at 
DREA[2,3j,  x  and  X  are  simply  proportional  to  one  another: 

*<*■•>  =  i  M 

where  x  runs  from  0  at  the  bow  to  20  at  the  stern  and  is  called  the  station  number.  The 
coordinate  a  runs  around  the  side  of  the  ship  and  has  values  in  the  range  [0,  l!.  The  keel  always 
has  the  value  a  =  1,  but  the  line  a  =  0  depends  upon  the  extent  of  the  hull  represented;  it  may 
be  the  sheerline,  the  design  waterline,  or  some  other  arbitrary  waterline.  If  the  hull  has  no 
knuckles,  a  is  proportional  to  the  arclength  around  each  station;  if  is  has  knuckles,  there  is  no 
simple  interpretation  of  the  coordinate  a.  The  coordinates  x  and  a  are  not  orthogonal.  Figure 
1  shows  the  two  different  coordinate  systems. 

The  functions  Y(x,  a)  and  Z(x,a)  have  the  form 

y  (x,  a)  =  By(x,  a)  ;  Z(x,  a)  =  Dz(x,  a)  (2.2) 
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Figure  1:  The  coordinate  systems 


where  B  is  the  half-breadth  of  the  hull  and  D  is  its  depth.  The  functions  y(x,  s)  and  z(x,  s ) 
are  complicated  and,  of  course,  are  different  for  every  hull  (see  Appendix  B).  They  are  given 
in  terms  of  a  large  number  of  spline  coefficients  (usually  several  thousand)  which  must  be 
determined  whenever  a  new  hull  is  to  be  represented.  These  coefficients  are  stored  in  a  data 
file.  Notice  that  the  dimensions  of  the  hull  can  be  changed  simply  by  changing  the  values  of  L, 
B,  and  D;  there  is  no  need  to  alter  the  functions  y(z,a )  and  z(x,  a).  This  is  the  simplest  form 
of  modification  allowed  by  MODHLL. 


3  Coordinate  Transformations 


Lackenby[4l  discusses  a  simple  method  by  which  hulls  can  be  modified  by  moving  the  hull 
sections.  The  hull  modification  is  done  by  performing  a  simple  coordinate  transformation  on 
the  coordinate  x.  This  amounts  to  a  shifting  of  each  section  forward  or  aft  by  an  amount 
that  is  dependent  on  x.  That  is,  if  z'  is  the  coordinate  on  the  original  hull  and  x(x')  on  the 
modified  hull,  then  station  x’  on  the  original  hull  is  shifted  to  station  x(x’)  on  the  modified  hull. 
The  forward  and  aft  perpendiculars  are  left  invariant  so  that  the  length  of  the  ship  between 
perpendiculars  remains  constant. 


Lackenby’s  method  involves  shifting  only  those  stations  forward  and  aft  of  any  parallel 
midbody.  Of  course,  shifting  the  stations  in  the  parallel  mid-body  would  cause  no  change  of 
the  ship’s  form.  The  coordinate  transformation  used  is  of  the  form 

x  =  a! xn  +  b'x'  for  x1  <  p^  (3.1) 

=  d*(20  -  x'f  +  e'(20  -  x ')  +  20  for  p'2  <  x'  <  20  (3.2) 

where  x  is  the  coordinate  describing  the  modified  hull,  s'  is  the  coordinate  describing  the 
original  hull,  p\  is  the  station  on  the  hull  at  which  the  parallel  mid-body  starts,  and  p2  is 
the  station  at  which  the  parallel  mid-body  ends.  If  no  parallel  mid-body  exists,  there  is  an 
implicit  assumption  that  the  tangents  to  the  hull  at  the  maximum  section  are  all  parallel  to  the 
centreplane,  otherwise  the  new  hull  would  not  be  smooth  at  the  mid-section.  This  method  is 
not  sufficiently  general  for  warship  hulls  which,  in  general,  will  not  satisfy  this  implicit  criterion. 
MODHLL  uses  a  generalization  of  Lackenby’s  method  which  relaxes  this  assumption. 

The  coordinate  transformation  used  is  similar  to  that  of  Lackenby,  but  extends  over  the 
mid-section  as  well  as  the  ends  of  the  hull.  It  is  then  no  longer  required  that  p\  and  p2  be  the 
start  and  end  of  parallel  mid-body:  they  may  be  any  two  stations  satisfying  0  <  p j  <  p2  <  20. 

If  the  old  hull  was  given  by  X'(x',a),  Y'( xf,a),  and  Z'{x',s),  the  functional  relations 
X(x,s),  Y(x,a),  Z{x,a )  for  the  modified  hull  are  given  by 


X(x,3)  =  X'(x\x),a) 

(3.3) 

Y(x,s)  =  r'(x'(x),s) 

(3.4) 

Z{x,  a)  =  Z\x’{x),a) 

(3.5) 

thus,  the  requirement  is  that  x'  be  given  in  term  of  x,  in  contrast  to  Lackenby’s  transformation. 

The  coordinate  transformation  used  is  represented  as  a  piecewise  cubic  polynomial  with 
breakpoints1  at  pi  and  pi ,  where  0  <  p\  <  p2  <  20;  pi  and  pi  are  otherwise  arbitrary.  If  p\  and 
P2  are  chosen  to  be  the  start  and  the  end  of  the  parallel  mid-body,  then  the  transformation 
assumes  the  same  form  as  that  of  Lackenby  except  for  the  reversal  of  roles  of  x  and  x'  and  the 
higher  order  of  the  polynomials  used.  Flow  calculations  which  use  the  hull  representation  some¬ 
times  require  that  X(x,  s),  Y(x,  a),  and  Z{x,a)  be  twice  differentiable.  This  in  turn  requires 
that  the  coordinate  transformation  have  two  continuous  derivatives  at  pi  and  p2 .  This  cannot 
be  done  unless  the  polynomial  pieces  are  of  at  least  fourth  order  (i.e.  cubic  polynomials).  Since 
there  are  three  polynomial  segments, 


x'(x)  =  a  +  bx  +  cx 2  +  dx3 

for  0  <  x  <  pi 

(3.6) 

=  e  +  fx  4-  gx2  +  hx3 

for  pi  <  x  <  p2 

(3.7) 

=  i  +  jx  +  kx*  +  lx3 

for  p2  <  x  <  20 

(3.8) 

there  are  12  coefficients  to  determine.  The  continuity  requirements  at  p\  and  p2  give  six  relations 
which  must  be  satisfied.  The  last  requirement  is  that  positions  of  the  bow  and  the  stern  do  not 
change;  thus  x  =  0  must  correspond  to  xf  =  0  and  x  =  20  to  x'  =  20.  This  gives  an  additional 

‘Breakpoints  are  the  places  at  which  the  piecewise  polynomial  segments  are  connected.  For  more  information 
on  numerical  calculations  using  piecewise  polynomials,  see  de  Boor[5] 
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two  conditions,  leaving  4  degrees  of  freedom  inherent  in  the  coordinate  transformation:  i.e. 
there  are  four  parameters  that  may  be  altered  when  specifying  hew  the  hull  is  to  change. 

The  transformation  is  defined  by  specifying  four  parameters  altered  in  one  of  the  following 
two  ways: 

1.  specifying  Xb,4>,p'\,  and  p'2,  or 

2.  specifying  4>p,4>a,Pi,  and  p^. 

<Pf  ,  4>a,  <?,  XB,p\  and  p2  are  defined  as  follows. 

1.  <j>r  is  the  prismatic  coefficient  of  the  forebody  measured  from  midships2. 

*'  =  To^io )/0  A(z)dz  (3'9) 

where  A(x)  is  the  curve  of  sectional  areas.  Note  that  A(x)  is  straightforward  to  calculate 
from  the  curve  of  sectional  areas  on  the  original  hull  since  A(x)  =  A'{x!{x)). 


2.  is  the  prismatic  coefficient  of  the  afterbody  measured  from  midships. 

3.  <t>  ^  (<I>a  +  4>f)/2,  the  total  prismatic  coefficient. 

4.  Xb  is  the  longitudinal  centre  of  buoyancy. 

1 


(3.10) 


XB  = 


20A(10)*;0 


rlQ 

/  xA{x)di 
Jo 


(3.11) 


5.  Pi  is  the  station  on  the  original  hull  which  correspond  to  p\.  This  may  be  used  to 
specified  the  start  of  a  parallel  mid-body.  Thus,  is  set  to  the  desired  start  of  the 
parallel  mid-body  on  the  the  new  hull,  and  pj  is  specified  to  be  the  start  of  the  parallel 
mid-body  on  the  original  hull. 

6.  pj  is  the  station  on  the  original  hull  which  correspond  to  pj.  This  may  be  used  to  specify 
the  end  of  a  parallel  mid- body. 

It  is  important  to  note  that  the  sectional  area  curve  for  a  modified  hull  is  always  calculated 
using  the  design  waterline  (see  Section  1),  even  if  the  draft  has  been  changed.  It  was  decided 
to  calculate  A(x)  in  this  way  for  three  reasons. 

3  Although  the  fore  and  aft  priematic  coefficient*  are  properly  meaeored  from  the  eection  of  maximum  area,  if 
they  are  measured  in  this  way,  change*  in  them  produce  little  change  in  the  form  of  hull*  with  near-parallel 
mid-body.  Hence,  for  the  purpoee*  of  this  memorandum,  it  i*  better  to  calculate  them  from  the  mid-section. 


1.  The  interpretation  of  station  0  and  20  as  the  forward  and  aft  perpendiculars  is  main¬ 
tained:  i.e.  they  are  the  stations  at  which  the  hull  surface,  the  centreplane  and  the  design 
waterline  intersect. 

2.  Changes  in  draft  will  not  affect  the  calculated  prismatic  coefficients.  To  maintain  constant 
prismatic  coefficients  when  the  draft  is  changed,  it  is  not  necessary  to  change  the  hull 
form,  as  would  be  the  case  if  the  sectional  areas  were  calculated  from  the  draft  waterline. 

3.  MODHLL  represents  the  draft  waterline  in  a  way  that  remains  invariant  if  the  trans¬ 
formation  x'(x)  is  changed  (see  Section  8).  Thus  the  draft  waterline  is  independent  of 
changes  in  hull  form.  However,  future  modifications  of  MODHLL  will  allow  calculation 
of  the  draft  waterline  for  different  states  of  ship’s  trim.  At  non-zero  trim  angles,  the  draft 
waterline  will  necessarily  become  dependent  on  changes  in  hull  form  as  well  as  the  trim 
angles  themselves.  If  the  curve  of  sectional  areas  were  to  be  calculated  using  the  draft 
waterline,  to  maintain  constant  prismatic  coefficients  while  changing  the  angle  of  trim 
would  require  a  change  in  hull  form.  Moreover,  calculation  of  the  transformation  x'(x) 
(see  Section  5)  would  be  greatly  complicated  by  the  fact  that  any  change  in  hull  form 
would  require  recalculation  of  the  draft  waterline. 


4  B-Splines 


Piecewise  polynomials  such  as  that  described  above  are  often  called  splines.  A  convenient 
way  of  ^presenting  a  spline  is  in  term  of  B-spline  basis  functions  (see  de  Boor[5]  for  an  ex¬ 
planation  of  B-splines).  This  allows  the  representation  of  the  transformation  in  the  following 
form: 

*'-*  =  £  **■(*)  (4-1) 

i=l 

where  Bi(x)  is  the  ith  B-spline  corresponding  to  the  knot  sequence  {0,0,0^1^2,20, 20, 20}1. 
By  choosing  only  three  identical  knots  at  each  endpoint  one  ensures  that  B,(0)  =  B,(20)  =  0 
and  hence  that  the  forward  and  aft  perpendiculars  remain  invariant  under  the  transformation. 
Figure  2  shows  the  four  B-splines  when  pi  =  8  and  P2  =  12.  The  spline  coefficients  £,  are 
adjusted  to  produce  the  transformation  desired.  If  the  change  in  the  hull  is  small,  then  each 
is  small:  i.e.  £,•  <  <  20. 

Changes  in  the  the  £,  cause  the  following  changes. 

‘A  knot  sequence  is  a  generalization  of  the  sequence  of  breakpoints;  each  element  of  the  sequence  is  called  a 
knot.  The  B-spline  basis  functions  depend  upon  the  positions  of  the  knots.  See  de  Boor(5j  for  an  exposition 
of  the  theory  of  splines. 
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If  each  of  the  specified  quantities  were  linear  in  the  the  calculation  of  the  &  would 
amount  to  the  solution  of  a  linear  system  of  equations  in  four  unknowns,  as  was  the  case  in 
Lackenby’s  method;  however,  due  to  the  inversion  of  the  roles  of  x  and  r1,  neither  the  prismatic 
coefficients  or  the  longitudinal  centre  of  buoyancy  are  linear  in  the  £,  since  they  appear  in  the 
argument  of  the  non-linear  function  A'(x')  used  to  evaluate  A(x).  The  method  for  calculating 
the  &  is  as  follows. 

Since  pj  and  pj  are  always  specified,  £i  and  £<  can  be  determined  in  terms  of  and  £3. 
Note  that  Bj(pj)  =  0  and  B+  (pi)  =  0.  Therefore, 


Pi  ~  Pi  =  CiBi(pi)  +  ?2-Bj(pi)  +  C3-Bs(pi) 
Pi  ~  Pi  =  +  &Bs(p2)  +  £4B4 (p2) 


£2  and  £3  are  found  by  a  simple  secant  search. 


/(x,y)  =  a 
y(x,y)  =  6 


be  two  simultaneous  equations  to  be  solved  for  x  and  y.  Expanding  to  first  order  around  x„ 
and  yn  one  gets 


/(**,  yn)  +  /*(*,  y)(x  -  In)  +  /y(x,  y)(y  -  y„)  «  a 
9{xn,  yn)  +  y*(x,  y)(x  -  xn)  -I-  gv{x,  y)(y  -  y„)  «  b 


where  the  superscripts  denote  partial  derivatives  by  x  or  y.  If  x„  and  y„  are  previous  approxi¬ 
mations  to  the  roots  of  equations  (5.3)  and  (5.4)  can  be  solved  to  yield  a  better  approximation: 


Xn+l  =  Xn  + 


gy(x,  y)(a  -  /(x,  y))  -  /y(x,  y)(b  -  g(x,  y)) 


yn+i  =  yn- 


/*(* ,  y)gv{x,  y)  -  /y(x,  y)y*(x,  y) 
9*{z,  y)(a  -  /(*,  y))  -  fx{x,  y){b  -  g{x,  y)) 
fz(x,y)g«(x,y)  -  fv{x,y)gx{x,y) 


The  derivatives  of  /  and  g  may  be  approximated  by  simple  finite  differences  of  previous  ap¬ 
proximations. 

The  following  algorithm  is  used  to  find  suitable  approximations  for  the  values  of  x  and  y. 


im 


Guess  xn,yn 

Calculate  /(xn,yn),y(xn,y„) 

Guess  xn+i 

Calculate  f(xn+i,yn),g(xn+i,yn) 

tz  /(*»+!> Vn)  ~  /{XntVn)  *  g(zn+l>yn)  —  g(zn>yn) 

/  = -  g  = - 

Xn+1  ~  Xn  Xn+l  “  Xn 

Guess  yn+i 

Calculate  /(x„+i ,  y„+i) ,  g (xn+i,  y„+i) 

Repeat 

n  =  n  +  1 

tu  /(®n*yn)  —  /(^-ni  yn-l)  u  g(*niyn)  ~~  g{xn,yn—i) 

J  =  -  g  =  - 

Vn  Vn—  1  Vn  S/n—  1 

_  _  ,  g>(g  -  /(xn,  yn))  -  /*(>  -  g(*n,  Vn)) 

n+1_  n+  f*9*-f*9* 

Calculate  /(x„+i,  yn),g(xn+i,  yn) 

^x  _  /(^n-H.Vn)  ~  f{*n,yn)  gX  _  gfcn-H,  y„)  ~  gpCn,  y„) 

Xn+1  —  Xn  zn+l  —  zn 


Vn+1  =  yn  + 


fx{b  -  g(xn+i,  y„))  -  gz(a  -  f[xn+l,yn)) 


fxgv  -  fyg‘ 

Calculate  /(xn+i ,  y„+i), g(*n+i,  Vn+l) 

Until  sufficient  accuracy  or  too  many  iterations 


When  4>F  and  4>a  are  specified,  the  independent  variables  of  both  functions  are  £2  and  £3. 
Thus,  substitute  the  following  values  into  the  preceding  algorithm: 


Xn 

-* 

&,n 

yn 

-* 

?3,n 

a 

— 

<t>F 

b 

- 

<t>A 

/(*»>  Vn) 

- 

foC&.n,  £s,n) 

g(xn,  yB) 

— 

4>A{^2,n,  & ,n) 

Whenever  4>f  and  4>a  are  to  be  evaluated  (i.e.  at  the  steps  requiring  the  calculation  of  /  and 
g),  £1  and  £4  must  first  be  determined  using  equations  (5.1)  and  (5.2).  <pp  and  4>a  are  then 
calculated  using  equations  (3.9)  and  (3.10). 

Similarly,  the  longitudinal  centre  of  buoyancy  and  the  total  prismatic  coefficient  also  de¬ 
pend  on  £2  and  £3.  Hence,  when  the  total  prismatic  coefficient  and  the  longitudinal  centre  of 
buoyancy  are  specified,  then  substitute  the  following  values  into  the  preceding  algorithm: 


a 


! 


f{xni Vn) 
}(ln, yn) 


6  Monotonicity  of  the  Transformation  | 

t 

The  function  describing  the  modified  hull  (x'(i)  =  x  +  £*=1  £,B,(x))  must  increase  with 
the  station  number,  z,  and  must  never  exceed  20  or  be  negative  because  the  functions 
*'(*', s),r(z',s)  and  are  not  defined  in  those  intervals.  Note  that  the  latter  two 

conditions  are  special  cases  of  the  first  one  because  the  knot  sequence  used  for  the  B-splines 

i  ensures  that  5,(0)  =  0  and  5,(20)  =  0.  < 

i 

j  Fritsch  and  Carlson[6]  have  shown  necessary  and  sufficient  conditions  for  a  cubic  inter-  ' 

J  polant  to  maintain  monotonicity.  Let  /(x)  be  the  cubic  interpolant,  let  (x,,  t  =  1, . . . ,  n}  be  its  j 

I  breakpoints  and 

A, 


<*i 

li 

sgn(a) 

Their  conditions  for  monotonic  data  sets  can  then  be  summarized  as  follows.  If  sgn(d,)  = 
sgn(dj+i)  =  sgn(Aj)  >  0,  then  /( z)  is  monotonic  on  [z<,Xj+j]  if  and  only  if  either 

1.  c»,-  +  7<  -  2  <  0;  or 

2.  a,-  +  "ft  —  2  >  0;  and  either 

(a)  2<Xi  -I-  7<  -  3  <  0;  or 


y 

(b)  a<  +  2~h  -  3  <  0;  or  \ 


:(*i) 


1  -  /, 

**  +  l  -  z« 

dA 

dxK 
Ai 

<U+i 

A, 

the  sign  of  a. 


-  XB 

Xsi&.n,  fs.n) 

Hb.n,  $3 ,n) 


This  test  is  used  to  check  the  monotonicity  of  the  transformation  between  the  modified  hull 
and  its  parent.  If  the  test  fails,  MODHLL  requires  either  that  the  modified  hull  be  discarded 
or  that  it  be  further  modified  until  i/( x )  is  monotonic. 


7  Calculation  of  Xb ,  and  4>a 


The  values  for  Xb,  4f,  and  4a  are  obtained  from  equations  (3.11),  (3.9),  and  (3.10),  each 
of  which  involves  the  curve  of  sectional  areas,  A(x).  This  function  is  itself  defined  in  terms  of 
a  line  integral  around  the  perimeter  of  section  x  as  follows. 


A(x)  =  -  f  YdZ  =  -  f\(x,»)^(x,s)ds 


Moreover,  the  curve  of  sectional  areas  of  a  modified  hull  is  given  by  A(x)  =  A'(x'(x))  where 
A'(x')  is  the  curve  of  sectional  areas  of  the  parent  hull.  Hence,  only  A'(x')  need  be  calculated. 
To  avoid  performing  this  integral  many  times,  A'(x')  is  calculated  once  for  21  values  of  x'  and 
the  resulting  values  are  splined.  The  resulting  spline  approximation  of  A'(x')  is  then  used  when 
calculating  Xb,  4f,  and  4a-  This  is  much  more  efficient  than  integrating  equation  (7.1)  many 
times  whenever  Xb  or  the  prismatic  coefficients  are  calculated. 

The  spline  of  the  curve  of  sectional  areas  is  calculated  using  the  algorithm  BSMTH  de¬ 
scribed  by  Hally  [7].  It  is  represented  in  terms  of  B-splines  as  follows. 

A'(x')=  't/3nBn(x')  (7.2) 

n=l 

where: 

•  j3n  is  the  nth  spline  coefficient, 

•  N  is  the  number  of  B-splines,  and 

•  Bn  is  the  nth  B-spline  of  fourth  order  corresponding  to  the  knot  sequence 
{0,0, 0,0, 2, 3, 4,..  .,16,17,18,20,20,20,20}. 

Each  of  the  above  integrals  is  calculated  numerically  by  sub-dividing  its  range  into  a 
number  of  smaller  intervals  and  using  four-point  Gaussian  quadrature  on  each  interval.  Since 
the  simple  linear  transformation  t  =  (2z  —  o  —  6) / (6  -  a)  will  translate  any  interval  'a,bj  into 
[-1,1]  provided  b  ^  a,  one  need  only  consider  integrals  over  the  interval  [-1,1]. 

Jm  f(x)dx  =  f_x  f  (— — a^-+~)  ^  dt  =  /  !  (7  3) 
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where 

,<0  =  /  (<^>L ±*±£ 

The  integral  of  y(t)  is  approximated  by 

l  » 

g(t)dt  as  ^  Ci9{ti)  (7.5) 

-1  .=1 

for  some  known  points  t,  and  coefficients  c<.  The  method  of  Gaussian  quadrature  chooses  these 
2n  parameters  optimally  so  that  all  polynomials  of  order  2n  or  less  are  calculated  exactly.  The 
U  for  which  this  occurs  are  the  roots  of  the  nth  order  Legendre  polynomial.  The  values  of  x, 
and  Ci  when  n  =  4  (used  by  MODHLL),  we  given  in  Table  1. 


(*-«) 

2 


(7.4) 


Roots,  ti 

Coefficients,  c< 

.8611363116 

.3399810436 

-.3399810436 

-.8611363116 

.3478548451 

.6521451549 

.6521451549 

.3478548451 

Table  1:  Roots  and  Coefficients  for  4- Point  Gaussian  Quadrature 


8  Representation  of  the  Draft  Waterline 


The  draft  waterline  of  the  hull  is  a  curve  *  =  w(x')  which  gives  the  s-value  of  the  draft 
waterline  as  a  function  of  the  station  number  on  the  parent  hull  x'.  The  function  tt^x')  is 
defined  implicitly  by 

Z'(x',w(x'))  =  Zw  (8.1) 

where  Zw  is  the  draft.  The  draft  waterline  on  a  modified  hull  is  given  by  s  =  u/(x'(x)).  Although 
the  HLLFLO  suite  of  programs  includes  procedures  for  solving  equation  (8.1)  for  w{x' )  (see 
Hally(8]),  they  are  not  very  efficient.  Hence,  a  spline  representation  of  the  curve  w(x')  is 
calculated  and  stored  in  the  hull  data  file,  so  that  the  s-value  of  the  draft  waterline  may  be 
obtained  quickly.  The  spline  is  calculated  using  the  algorithm  BSMTH  described  by  Hally  [7], 
It  is  represented  in  terms  of  B-splines  as  follows. 

N 

w(x')  =  £  OnBnA*')  (8.2) 

n=  1 


where: 


y  t 


•  /3n  is  the  nth  spline  coefficient, 

•  N  is  the  number  of  B-splines, 

•  A:  is  the  order  of  the  spline  (usually  4),  and 

•  Bnj t  is  the  ntfl  spline  corresponding  to  a  set  of  knots  in  which  there  are  k  knots  at  each 
end-point  of  the  draft  waterline  and  the  remaining  (N  +  k)  —  2k  knots  are  equally  spaced. 

If  the  depth  of  the  hull,  D,  is  changed  but  Z w  remains  constant,  the  draft  waterline 
function  tv(x')  also  changes.  To  avoid  the  recalculation  of  w(x'),  whenever  D  is  changed  by  a 
factor  a  (i.e.  D  — ►  aD)  Zw  is  also  scaled  by  a  (i.e.  Zw  — *•  aZw).  In  this  case,  tu(x')  remains 
invariant  so  that  the  waterline  curve  need  not  be  recalculated  or  resplined. 


9  Concluding  Remarks 


In  order  to  see  how  the  flow  and  the  resistance  around  a  ship  change,  it  is  useful  and 
convenient  to  generate  a  slightly  altered  hull  from  one  whose  computer  representation  is  already 
defined.  The  MODHLL  program  uses  a  generalization  of  Lackenby’s  method  on  which  the  new 
hull  is  created  by  specifying  the  new  characteristics  of  the  hull  such  as  the  prismatic  coefficients 
and  centre  of  buoyancy.  Furthermore',  it  is  also  possible  to  change  the  hull’s  dimensions  and  its 
draft  without  affecting  the  basic  characteristics  of  the  hull.  The  modified  hull  is  defined  from 
the  parent  hull  by  means  of  a  simple  spline  transformation  relating  the  station  numbers  x  and 
x'  on  the  respective  hulls. 

For  the  sake  of  efficiency,  the  program  uses  a  spline  representation  to  define  the  sectional 
area  curve  all  along  the  hull.  This  curve  gives  the  sectional  area  at  any  station  on  the  hull  and 
is  used  for  calculating  the  prismatic  coefficients  and  the  centre  of  bouyancy  of  both  modified 
and  parent  hulls. 
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This  appendix  u  a  user’s  guide  for  the  program  MODHLL.  MODHLL  has  been  designed 
to  be  integrated  within  the  HLLFLO  suite  of  programs  and  therefore  many  implementation 
decisions  have  been  made  to  provide  consistency  with  the  other  programs  in  HLLFLO.  The 
first  five  chapters  of  the  HLLFLO  User’s  Guidejl]  provide  overall  information  about  HLLFLO. 

MODHLL  itself  is  written  in  ANSI  standard  FORTRAN  77  but  calls  HLLFLO  proce¬ 
dures  which  have  minor  non-ANSI  features:  see  the  HLLFLO  User’s  Guide,  Chapter  5  for 
details.  MODHLL  calls  procedures  from  the  spline  libraries  HLLYSP[7]  and  BSPLIN[5],  from 
the  GETWRD  Package[9],  and  from  the  PLOTlOflO]  or  IGL[11]  graphics  libraries.  In  addition, 
procedures  are  called  from  the  HLLFLO  libraries  CMDSUB,  HCRSUB,  HLLSUB,  MSGSUB, 
PLTSUB,  and  SMVSUB. 

MODHLL  is  designed  to  be  run  on  a  Tektronix  graphics  terminal  or  equivalent.  Its 
graphics  are  generated  by  the  HLLFLO  PLTSUB  library  which  in  turn  call  procedures  from 
either  the  PLOTlO  library  or  the  IGL  library.  See  the  HLLFLO  User’s  Guide,  Chapter  4  for 
details. 

On  the  DEC- 20/60  at  DREA,  a  simple  method  for  running  MODHLL  is  available.  A  new 
EXEC  command  may  be  defined  by  typing 

(BDECLARE  PCL  PS :  <HALLY>HLLFLQ 

at  the  EXEC  level  of  command.  This  uses  the  PCL  extension  to  the  TO  PS-20  EXEC  to  define 
the  new  command  HLLFLO.  By  giving  the  command 

0HLLFL0  MODHLL  plot-library 

MODHLL  is  automatically  compiled,  loaded,  and  run.  If  the  file  MODHLL.EXE  already  exists 
on  the  directory  to  which  you  are  connected,  it  will  be  run.  Otherwise  a  search  is  made  for  .REL 
versions  of  all  necessary  library  files.  If  an  .REL  version  of  a  file  cannot  be  found,  the  source 
file  on  PS:<HALLY>  is  compiled.  Once  .REL  versions  of  all  necessary  files  have  been  created, 
they  are  linked  to  create  the  file  MODHLL.EXE  which  is  then  run.  New  .REL  and  .EXE  files 
remain  on  the  connected  directory  for  subsequent  use.  The  sub-command  plot-library  is  one  of 
PLOTlO  or  IGL  according  to  the  plotting  library  you  wish  to  use. 

MODHLL  uses  a  command  language  to  allow  you  to  define  how  hulls  are  to  be  modified: 
see  Chapter  3  of  the  HLLFLO  User’s  Guide[l]  for  details.  The  prompt 

command«> 

indicates  that  MODHLL  is  awaiting  command.  The  allowed  commands  and  their  functions  are 
described  in  the  following  sections.  The  command 


command“>EXIT 


is  used  to  stop  execution  of  the  program. 


A.l  Defining  the  Terminal  Characteristics 

Upon  entering  MODHLL  the  message 
Entering  MODHLL  . . . 

terminal  type?  *> 

is  written  on  the  terminal  screen.  This  is  the  prompt  to  identify  the  characteristics  of  the 
terminal  which  you  are  using.  See  the  HLLFLO  User’s  Guide,  Chapter  4  for  an  explanation  of 
the  possible  responses.  The  terminal  characteristics  may  be  changed  by  giving  the  command 

c  ommand*>TERMINAL,‘ 

followed  by  the  characteristics  of  the  terminal.  Typically  you  will  plot  a  graph  in  MODHLL 
then  use  this  command  to  redefine  the  plotter  to  be  the  laser-printer  rather  than  the  terminal 
so  that  a  hard  copy  of  the  graph  can  be  obtained. 


A. 2  Loading  a  Parent  Hull  File 

At  all  times  MODHLL  saves  the  characteristics  of  two  hulls:  the  parent  and  a  modified 
form  of  it.  Before  modifications  of  a  hull  can  be  made,  a  parent  hull  must  first  be  loaded  into 
the  program  memory.  A  prompt  for  the  name  of  the  file  containing  the  parent  hull  data  occurs 
when  you  enter  MODHLL,  immediately  after  the  prompts  for  the  terminal  characteristics. 
Subsequently  a  new  parent  hull  can  be  loaded  by  giving  the  command 

command=>LOAD 

after  which  you  will  be  prompted  the  name  of  the  parent  hull  data  file.  The  file  name  may  be 
at  most  30  characters  long  including  directory  specifications  and  extensions.  If  no  extension 
is  given,  the  default  on  the  DREA  DEC-20/60  is  .DAT.  The  input  file  must  be  of  the  same 
format  as  an  output  file  from  the  HLLFLO  program  SMHULL  (see  Appendix  B)  or  an  output 
file  from  MODHLL  (see  Section  A. 3).  In  the  latter  case,  the  name  of  the  file  containing  the 
parent  hull  data  is  read  from  the  file  containing  the  modified  hull  data. 

When  a  parent  hull  is  loaded,  the  modified  hull  is  set  to  be  the  same  as  the  parent.  This 
causes  all  characteristics  of  the  previous  modified  hull  to  be  lost  from  the  program  memory.  If 
you  have  created  a  modified  hull  and  then  give  the  LOAD  command  before  saving  it  in  a  file  (see 
Section  A. 3),  you  are  asked  whether  you  wish  to  save  the  previous  modified  hull  before  loading 
the  new  parent.  This  prevents  loss  of  time  if  you  carelessly  forget  to  save  the  modified  hull  you 
have  just  created. 


STATION 


Figure  A.l:  Spline  of  Normalized  Curve  of  Sectional  Areas 

Whenever  a  parent  hull  is  loaded,  a  spline  representation  of  its  curve  of  sectional  areas  is 
calculated  (see  Section  7).  A  message  is  written  on  the  screen  when  the  calculation  begins  and 
ends.  You  are  then  asked  whether  you  would  like  to  examine  a  display  of  the  spline  curve.  It 
is  usually  wise  to  do  so  since  splines  will  sometimes  contain  unwanted  wiggles  that  could  cause 
problems  when  calculating  the  prismatic  coefficients.  If  you  answer  YES,  a  display  such  as  that 
in  Figure  A.l  is  plotted.  Notice  that  it  is  actually  the  curve  A'{x')/ BZj,  that  is  splined,  where 
B  is  the  half-breadth  and  Z*  is  the  design  waterline  of  the  parent  hull. 

When  the  display  has  been  plotted  the  prompt 

sub-command*> 

appears  in  the  top  left  corner  of  the  screen.  If  the  spline  curve  is  not  to  your  liking  it  may  be 
altered  as  described  in  Section  A. 8.  To  accept  the  spline  curve  use  the  sub-command  CONTINUE. 
The  screen  will  be  cleared  and  the  characeristics  of  the  parent  hull  will  be  displayed.  MODHLL 
then  awaits  further  commands. 

A. 3  Saving  a  Modified  Hull 

A  modified  hull  may  be  saved  by  giving  the  command 
command«>SAVE 

You  are  prompted  for  a  title  for  the  hull  (for  use  on  graphic  displays  generated  by  HLLFLO 
graphics  programs),  and  for  the  name  of  the  file  in  which  to  store  the  data.  The  former  is  a 


string  of  at  most  60  characters,  the  latter  a  string  of  at  most  30  characters. 

You  cannot  save  a  hull  whose  x-transformation  has  not  yet  been  calculated  (since  this 
transformation  is  integral  to  the  definition  of  the  modified  hull),  nor  can  you  save  a  hull  whose 
x-transformation  is  not  monotonic. 

The  file  created,  which  contains  all  relevant  information  defining  a  modified  hull,  can  be 
read  subsequently  by  other  programs  (for  example,  a  program  calculating  the  flow  around  this 
modified  hull). 

The  modified  hull  data  file  is  written  in  the  following  format. 


Variable 

Format 

‘MODIFIED  HULL’ 

(A) 

HLLFIL 

(A) 

TITLE 

(A) 

SLNGTH,  BDTH,  DPTH 

(3G15.7) 

PI,  P2 

(2G15.7) 

XI(I),  1=1,..  .,4 

(4G15.7) 

If  ZWTR  ^  DPTH-. 

ZWTR 

(G15.7) 

KWTR,  NWTR 

(215) 

TWTR(J),  J=1,...,KWTR+NWTR 

(5G15.7) 

BWTR(J),  J=1 . NWTR 

(5G15.7) 

where: 

•  ‘MODIFIED  HULL’  is  a  character  string  flag  used  to  indicate  that  the  file  contains  data 
describing  a  modified  hull  rather  than  a  parent, 

•  HLLFIL  is  the  name  of  the  file  containing  the  parent  hull  data, 

•  SLNGTH,  BDTH  and  DPTH  are  respectively  the  length,  the  half  breadth  and  the  depth 
of  the  ship, 

•  Pi  is  the  station  number  p j, 

•  P2  is  the  station  number  pi, 

•  XI  is  an  array  containing  the  values  of  &,  i  =  1, . . . ,  4, 

•  ZWTR  is  the  draft  of  the  ship, 

•  KWTR  is  the  order  of  the  spline  used  to  define  the  position  of  the  draft  waterline, 

•  NWTR  is  the  number  of  B-splines  used  to  define  the  position  of  the  draft  waterline. 

•  TWTR  is  the  knot  sequence  used  to  define  the  position  of  the  draft  waterline,  and 

•  BWTR  are  the  spline  coefficients  used  to  define  the  position  of  the  draft  waterline. 


A. 4  Modifying  the  Hull 


As  described  in  Section  A. 2,  MODHLL  keeps  the  characteristics  of  a  parent  hull  and  a 
modified  hull  in  its  memory.  The  modified  hull  starts  out  as  an  exact  replica  of  the  parent  but 
can  be  modified  by  the  commands  described  in  this  section.  The  parent  hull  remains  invariant 
until  a  new  parent  is  loaded  using  the  LOAD  command.  The  command  to  change  the  modified 
hull  is 

command*>CHANGE  parameter  value 

where  parameter  is  one  of  AFT-PRISMATIC-COEF,  DEPTH,  DRAFT,  FORWARD-PRISMATIC-COEF, 
HALF-BREADTH,  LCB,  LENGTH,  PI,  PI  \  P2,  P2*,  or  TOTAL-PRISMATIC-COEF  and  value  is  the  new 
value  far  that  parameter.  The  default  value  for  p\  and  p[  is  8.0.  The  default  value  for  p2  and 
p'2  is  12.0.  The  default  values  for  the  remaining  parameters  are  the  corresponding  values  on 
the  parent  hull. 

A.4.1  Changes  in  the  Draft 

When  the  draft  of  the  hull  is  changed,  a  new  draft  waterline  curve  w(x')  must  be  calculated. 
As  described  in  Section  8,  w(x')  is  approximated  by  a  spline.  Before  the  spline  can  be  calculated, 
its  order  (the  order  of  the  polynomial  segments  of  which  it  is  comprised)  and  the  number  of 
B-splines  used  to  represent  it  must  be  specified  (see  de  Boor[5]  for  further  information  about 
splines).  Reasonable  values  for  these  variables  are  4  and  20  respectively.  The  values  are  input 
as  sub-commands  to  the  command  CHANGE  DRAFT.  The  full  command  has  the  form 

c  o  mm  and”  >  CHANGE  DRAFT  value-of-draft  order-of-spline  no-of-B-splines 

so  that  a  typical  command  might  be 

c ommand” > CHAN GE  DRAFT  15.5  4  20 

After  the  command  has  been  given,  MODHLL  calculates  the  spline  and  then  displays  it  so  that 
you  can  examine  it.  When  the  display  has  been  plotted  the  prompt 

sub-command=*> 

appears  in  the  top  left  corner  of  the  screen.  If  the  spline  curve  is  not  to  your  liking  it  may  be 
altered  as  described  in  Section  A.8.  To  accept  the  spline  curve  use  the  sub-command  CONTINUE. 
The  screen  will  be  cleared  and  MODHLL  then  awaits  further  commands. 

A. 4. 2  Changes  which  Affect  x'(x) 


A  change  in  the  value  of  the  length,  half-breadth,  depth,  or  draft  of  the  hull  in  no  way 
changes  the  transformation  z'(x).  If  only  these  parameters  are  changed  the  values  of  £,  which 
define  this  transformation  need  not  be  recalculated;  however,  a  change  in  any  of  the  other 
parameters  requires  that  the  be  recalculated  before  the  modified  hull  can  be  saved  or  its 
curve  of  sectional  areas  can  be  plotted.  The  command  to  calculate  the  transformation  is 


command->CALCULATE-X-TRANSFORMATIOH 


Note  that  you  can  make  as  many  changes  to  the  hull  as  you  wish  before  calculating  the  £,. 

As  described  in  Section  3,  the  x-transformation  can  be  calculated  using  the  values  of  either 
Xb  and  4>  or  of  <Pf  and  <p, If  only  the  values  of  p\,  p'1(  pj,  and  p2  are  changed,  the  default  is  to 
use  Xg  and  <t>.  Otherwise,  MODHLL  will  determine  which  parameters  have  been  changed  and 
calculate  the  i-transformation  appropriately.  For  example,  if  only  <f>p  has  been  changed,  <j>p 
and  <j>A  will  be  used.  In  the  case  of  an  ambiguity  (e.g.  if  both  Xb  and  <pA  have  been  changed), 
you  will  be  prompted  to  clarify  which  set  of  parameters  you  wish  to  be  used. 

Whenever  an  i-transformation  is  calculated  successfully,  the  parameters  of  the  modi¬ 
fied  hull  described  by  that  transformation  become  the  new  defaults.  For  example,  suppose 
the  current  default  value  for  Xb  is  10.4  and  of  4>p  is  0.61.  You  change  <pF  to  0.59  and  the 
i-transformation  is  calculated  successfully.  MODHLL  recalculates  Xb  using  the  new  trans¬ 
formation  obtaining  a  value  of  10.4782.  This  becomes  the  new  default  value  for  Xb  and  the 
default  for  4>p  is  0.59. 

When  calculating  the  MODHLL  prints  on  the  terminal  a  record  of  the  number  of  iter¬ 
ations  performed  so  far.  If  ten  iterations  ewe  performed  without  convergence,  the  calculation  is 
ended.  MODHLL  stores  the  values  of  the  which  gave  the  most  accurate  fit  to  the  parameters 
you  specified.  You  are  given  the  option  of  retaining  the  hull  as  defined  by  these  “best”  ^-values 
as  the  new  modified  hull.  Remember  that  if  you  do  so,  the  default  values  for  all  parameters 
will  be  set  to  the  values  for  this  hull.  If  you  do  not  accept  the  “best”  f-values,  you  must  make 
further  changes  and  recalculate  the  transformation. 

Sometimes  when  an  attempt  to  calculate  the  i-transformation  fails,  you  will  wish  to  restore 
all  defaults  to  the  values  for  the  parent  hull.  This  may  be  done  by  giving  the  command 

c  ommand">RESET 

As  noted  in  Section  4,  the  values  of  p\,  p\,  pi,  and  p2  affect  the  values  of  the  prismatic 
coefficients  for  the  modified  hull  as  compared  with  their  values  for  the  parent  hull.  The  most 
common  cause  of  failure  in  the  calculation  of  the  x-transformation  is  to  specify  values  of  the 
p’s  which  tend  to  increase  a  prismatic  coefficient,  while  requiring  that  the  prismatic  coefficient 
actually  decrease,  or  vice  versa.  For  example,  setting  pi  <  p\  while  requiring  that  4>f  decrease 
will  likely  cause  trouble.  Before  starting  the  calculation  of  the  i- transformation,  MODHLL 
checks  the  values  of  the  p’s  and  of  the  required  changes  in  the  primatic  coefficients.  If  it 
suspects  there  may  be  difficulties,  you  will  be  notified  of  the  problem  and  asked  whether  you 
wish  to  continue  with  the  calculation  anyway.  The  following  message  is  typical. 

command*>CALC 

P2-P1  »  4.000  and  P2’-P1‘  *  3.000  tending  to  increase  the  total 
prismatic  coefficient.  But  you  wish  to  decrease  the  total 
prismatic  coefficient.  The  iteration  in  the  calculation  of 
the  x*transformation  is  likely  to  diverge. 

Do  you  wish  to  proceed  with  the  calculation  (Y  or  H)? 


19 


Table  2  summarizes  the  changes  suggested  in  the  p’s  to  avoid  difficulties  when  changing  the 
hull  parameters. 


Parameter 

Modification 

Suggested 
Change  in  p’s 

Xg  >  X'R 

Pi  +  Pi>  pi;  +  p\ 

<t>  >  <t>' 

Pi  ~  Pi  >  Pa  ~  Pi 

<t>p  >  4>'F 

Pi  <  Pi 

<t>A  >  4>'a 

Pi  >  p'2 

Table  2:  Summary  of  Suggested  Changes  in  pi  and  p2 


A. 5  Printing  Hull  Characteristics 


After  changing  a  few  hull  parameters,  it  is  easy  to  forget  the  value  of  a  specific  parameter 
on  the  modified  or  parent  hull.  The  values  of  all  parameters  can  be  printed  on  the  terminal 
using  the  commands 

command«>PRINT-MODIFIED-HULL 

and 


c  ommand*  >PRINT - PARENT - HULL 

The  following  is  an  example  of  output  generated  by  PRINT-MOD IFIED-HULL. 


The  characteristics  of  the  current  hull  are: 


PI’  -  8.50  PI  *  8.00 

P2 '  *  11.00  P2  ■  12.00 

Longitudinal  centre  of  buoyancy 
Total  prismatic  coefficient 
Forward  prismatic  coefficient 
Aft  prismatic  coefficient 
Length  between  perpendiculars 
Half  breadth 
Depth 
Draft 


10.3732 
0.5917 
unknown 
unknown 
112.000  M 
7.260  M 
11.500  M 
8.125  M 


The  unknown  values  of  4>p  and  4>a  occur  because  either  Xg  and  or  <j>  has  been  changed  but  the  x- 
transformation  has  not  yet  been  recalculated,  so  that  the  values  for  4>f  and  4>a  which  correspond 
to  the  given  values  of  Xg  and  ^  have  not  yet  been  determined.  If  the  r-transformation  had 
already  been  determined,  the  values  for  fo  and  would  be  included. 
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A. 6  Plotting  the  Curves  of  Sectional  Areas 


A  display  of  the  curve  of  sectional  areas  of  both  the  modified  hull  and  the  parent  hull  can 
be  generated  by  giving  the  command 

c  ommand«>PLOT-SECTIONAL- AREAS 

An  example  of  the  display  is  shown  in  Figure  A.2. 

Note  that,  even  if  the  draft  has  been  modified,  the  area  of  each  section  is  defined  to  be  the 
area  below  the  design  waterline:  the  waterline  of  the  parent  hull  perhaps  scaled  by  a  change 
in  depth  (see  Sections  1  and  8).  For  example,  suppose  the  depth  of  the  parent  hull  is  equal  to 
20  ft  and  its  design  waterline  is  15  ft  above  the  baseline.  If  the  depth  of  the  ship  is  modified 
to  a  value  of  25  ft,  then  the  design  waterline  for  the  modified  hull  is  18.75  ft  from  the  baseline. 
Even  if  the  draft  of  the  modified  hull  is  subsequently  changed  to  17  ft,  the  sectional  areas  are 
still  calculated  using  the  value  18.75  ft. 


CURVE  OF  SECTIONAL  AREA  FOR 
PARENT  AND  MODIFIED  HULL 


A. 7  Plotting  x'(x ) 

A  display  of  the  currently  defined  transformation  function  x'(x)  can  be  generated  by  giving 
the  command 

command*>PL0T-X-TRANSF0RMATI0N 
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An  example  of  the  display  is  shown  in  Figure  A.  3.  The  straight  line  on  the  display  indicates 
the  transformation  which  defines  an  unmodified  hull  (z'(z)  =  z). 


TRANSFORMATION  CURVE  FOR 
MODIFIED  HULL 


Figure  A.3:  Transformation  curve  for  a  modified  hull 


A.8  Modification  of  Splines 

When  a  spline  of  either  the  sectional  area  curve  or  the  design  waterline  of  the  parent  hull 
has  been  calculated  and  plotted,  MODHLL  enters  a  sub-command  mode  which  allows  you  to 
alter  the  spline  if  it  is  not  to  your  liking.  Figure  A.l  is  an  example  of  the  appearance  of  the 
terminal  display  at  this  point.  The  positions  of  the  knots  used  for  the  spline  are  shown  by  the 
crosses  near  the  bottom  of  the  display. 

The  prompt 
aub-command*> 

which  appears  at  the  top  left  of  the  screen  indicates  that  MODHLL  is  awaiting  your  command. 
If  the  spline  is  satisfactory,  use  the  sub-command 

sub-command*>CONTINUE 

to  continue;  however,  if  you  wish  to  make  changes  to  the  spline,  use  the  command 
8ub-command*>CHANGE  parameter 

one  or  more  times  to  alter  the  spline  parameters.  The  sub-command 
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sub-command*>RESPLINE 


will  then  recalculate  and  display  the  modified  spline. 

The  following  sections  describe  the  spline  parameters  which  may  be  changed  in  order  to 
modify  the  spline  curve. 


A. 8.1  Changing  the  Amount  of  Smoothing 


The  splines  calculated  by  MODHLL  are  smoothing  splines  and  hence,  the  desired  amount 
of  smoothing  must  be  specified.  (The  spline  is  calculated  by  the  HLLYSP  library  subroutine 
BSMTH  whose  algorithm  is  discussed  by  Hally (7].)  This  is  done  by  establishing  an  appropriate 
X2  for  the  spline  defined  by 


2  V-*  (/(*»)  ~  j/n)2 

X  .2 


(A.l) 


where  (zn, yn), n  =  1 ,N  are  the  data  points,  f{x)  is  the  function  which  defines  the  spline, 
and  en  is  an  error  associated  with  the  nth  data  point.  The  accuracy  of  the  spline  is  defined  to  be 
\Zx2/N-  If  all  the  e„  are  1.0,  then  the  accuracy  of  the  spline  is  the  same  as  the  mean  distance 
of  the  data  points  from  the  spline  curve.  Conversely,  if  the  accuracy  is  1.0,  the  distance  of  the 
spline  from  each  data  point  will  be  on  average  equal  to  its  error.  The  accuracy  of  the  spline 
may  be  changed  using  the  command 


conun  and->CHANGE  AMOUNT  -  OF-  SMOOTHING  method 


where  method  is  one  of  FOR-ALL-SPLINES-  or  CALCULATED-FROM-DATA. 

The  command  FOR-ALL-SPLINES-  requires  an  argument:  a  number  equal  to  the  accuracy 
for  each  of  the  splines. 

The  command  CALCULATED-FROM-DATA  indicates  that  the  accuracy  of  each  spline  is  to  be 
calculated  internally  by  the  function  PRERR  in  the  HLLYSP  library [7].  This  is  the  default. 
Once  changed  the  default  for  the  amount  of  smoothing  remains  at  its  new  value  until  changed 
again.  If  the  accuracy  is  calculated  from  the  data  by  PRERR,  all  subsequent  splines  will  also 
call  PRERR  until  the  sub-command  CHANGE  AMOUNT-OF-SMOOTHING  is  invoked  again. 


A. 8. 2  Changing  the  Data  Point  Errors 

The  data  point  errors,  e„,  may  also  be  changed.  The  smaller  the  error,  the  closer  the 
spline  will  try  to  come  to  its  corresponding  data  point.  The  appropriate  command  is 

sub-command*>CHANGE  ERRORS  method 

where  method  is  one  of  ALL-EQUAL  or  INTERACTIVE-SPECIFICATION.  The  command  ALL-EQUAL 
causes  the  value  of  all  errors  to  be  set  to  1.0. 

The  command  INTERACTIVE-SPECIFICATION  allows  the  changing  of  the  errors  of  individ¬ 
ual  data  points  via  the  cross-hair  sub-command  mode  (see  Section  3.1  of  the  HLLFLO  User’s 
Guidefl]).  Upon  giving  the  command,  the  data  points  of  the  current  spline  are  replotted  and 


the  cross-hair  appears  as  a  prompt  for  a  single  letter  command.  The  available  commands  are 
the  following. 

1.  E  or  e:  Exits  from  the  cross-hair  sub-command  mode  to  the  sub-command  mode 

2.  N  or  n:  A  new  value  for  the  error  of  the  data  point  closest  to  the  cross-hair  may  be 
entered.  A  line  is  drawn  from  the  data  point  identified  to  the  value  typed  in.  The  line 
allows  you  to  remember  which  error  value  is  associated  with  which  data  point. 

3.  O  or  o:  The  old  value  of  the  error  of  the  data  point  closest  to  the  cross-hair  is  printed  on 
the  screen  somewhat  above  or  below  the  cross-hair  position. 

4.  P  or  p:  Plots  the  data  points  again.  After  several  errors  have  been  changed,  the  screen 
can  get  cluttered  with  old  and  new  error  values.  This  command  allows  the  screen  to  be 
cleared  of  all  these  numbers. 

5.  ?:  Causes  a  list  of  possible  cross-hair  sub-commands  to  be  written  above  the  plot. 

The  errors  once  changed  retain  their  altered  values  until  changed  again. 

A. 8. 3  Changing  the  Stiffness  Weights 

The  stiffness  weights  are  parameters  which  can  be  used  to  change  the  “stiffness”  of  the 
spline  in  the  interval  between  any  two  successive  knots.  By  increasing  the  stiffness  weight 
between  two  knots,  the  spline  can  be  made  straighter  there.  By  decreasing  the  stiffness  weight 
in  a  given  knot  interval,  the  spline  can  be  allowed  to  bend  more  easily  in  that  interval.  Thus, 
it  is  desirable  to  have  large  stiffness  weights  where  the  spline  should  vary  very  little,  and  small 
stiffness  weights  where  the  spline  is  varying  rapidly.  To  change  the  stiffness  weights  give  the 
command 

command*>CHANGE  STIFFNESS-WEIGHTS  method 

where  method  is  one  of  ALL-EQUAL,  INTERACTIVE-SPECIFICATION  or  CALCULATED-FROM-DATA. 

The  command  ALL-EQUAL  sets  all  the  stiffness  weights  to  1.  Since  only  the  relative  values 
of  the  stiffness  weights  are  important,  there  is  no  need  to  allow  the  setting  of  all  stiffness  weights 
to  a  value  other  than  1. 

The  command  INTERACTIVE-SPECIFICATION  allows  the  changing  of  individual  stiffness 
weights  via  the  cross-hair  sub-command  mode  (see  Section  3.1  of  the  HLLFLO  User’s  Guide|l]). 
Upon  giving  the  command,  the  data  points  of  the  current  spline  are  replotted  and  the  cross-hair 
appears  as  a  prompt  for  a  single  letter  command.  The  knots  for  the  spline  are  indicated  by  the 
row  of  crosses  at  the  bottom  of  the  plot.  There  is  a  stiffness  weight  associated  with  the  interval 
between  each  pair  of  adjacent  knots.  The  available  commands  are 

•  E  or  e:  Exits  from  the  cross-hair  sub-command  mode  to  the  sub-command  mode. 
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•  N  or  n:  A  new  value  for  the  stiffness  weight  of  the  knot  interval  identified  by  the  vertical 
cross-hair  may  be  entered.  A  line  is  drawn  from  the  centre  of  the  knot  interval  to  the 
position  of  the  cross-hair.  The  line  allows  you  to  remember  which  value  is  associated  with 
which  knot  interval. 

•  O  or  o:  The  old  value  of  the  stiffness  weight  of  the  knot  interval  identified  by  the  vertical 
cross-hair  is  printed  at  the  position  of  the  cross-hair. 

•  P  or  p:  Plots  the  data  points  again.  After  several  stiffness  weights  have  been  changed, 
the  screen  can  get  cluttered  with  old  and  new  stiffness  weight  values.  This  command 
allows  the  screen  to  be  cleared  of  all  these  numbers. 

•  ?:  Causes  a  list  of  possible  cross-hair  sub-commands  to  be  written  above  the  plot. 

The  command  CALCULATED -FROM-DATA  indicates  that  the  stiffness  weights  for  each  spline 
are  to  be  calculated  by  the  subroutine  WTIBEG  in  the  HLLYSP  library [7],  This  is  the  default. 

New  values  for  the  stiffness  weights  remain  as  defaults  until  changed  again.  If  the  stiffness 
weights  are  calculated  from  the  data  by  WTIBEG,  all  subsequent  splines  will  also  call  WTIBEG 
until  the  command  CHANGE  STIFFNESS -WEIGHTS  is  invoked  again. 

A. 9  Spline  Trouble-shooting 

MODHLL  calculates  the  spline  coefficients  for  each  spline  by  evaluating  and  inverting  a 
matrix.  Occasionally,  the  matrix  will  be  nearly  singular  and  the  resulting  spline  will  be  wild. 
MODHLL  is  almost  always  able  to  detect  the  problem.  If  the  spline  is  not  plotted,  the  message 

??Error  while  calculating  spline. 

Probably  numerical  difficulties. 

is  written  on  the  terminal  screen.  If  the  spline  is  plotted,  no  message  is  written,  but  it  it  will 
be  evident  from  the  plot  that  there  is  a  problem.  Often  the  values  of  the  spline  will  extend  to 
huge  numbers:  e.g.  1010  or  1020.  It  is  usually  possible  to  fix  the  spline  by  changing  the  amount 
of  smoothing.  In  the  typical  case,  the  required  accuracy  of  the  spline  is  too  low:  try  increasing 
it  by  a  factor  of  10  or  100  and  see  if  that  works.  To  determine  the  least  amount  of  smoothing 
possible,  you  can  later  decrease  it  again  progressively  until  the  numerical  problems  reappear. 
Numerical  difficulties  can  also  be  caused  if  the  values  of  the  errors  or  of  the  stiffness  weights 
are  set  too  high  or  too  low.  The  most  stable  situation  is  if  all  errors  and  all  stiffness  weights 
are  set  equal  to  1.0. 


Appendix  B  Modifications  to  the 
HLLFLO  Hull  Representation 


In  order  to  facilitate  the  modification  of  hull  dimensions  by  MODHLL,  the  HLLFLO  hull  B- 
spline  representation  has  been  modified  slightly.  In  the  previous  representation  (see  Hally  [2,3] ), 
the  cartesian  coordinates  ( X ,  Y,  Z)  were  given  as  function  of  the  non-cartesian  coordinates  (x,  s) 
by 
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where  B^ki(x),n  =  I ,  ...,iVx,  are  B-spliaes  of  order  fcx  corresponding  to  a  knot  sequence 
tL*\  n  =  1 ,  ■  ■  ■ ,  Nx  +  kx,  and  Bf'kt(a),j  =  1 ,N,  are  B-splines  of  order  ks  corresponding  to 

the  knot  sequence  t^,  n  =  1, . . . ,  N,  +  k,.  In  the  modified  scheme,  the  spline  coefficients  are 
non-dimensionalized  using  the  hull  half-breadth,  B,  and  depth,  D 
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The  virtue  of  the  new  scheme  is  that  none  of  the  spline  coefficients  need  be  changed  when  B 
and  D  are  changed.  In  the  old  scheme,  if  B  were  doubled,  then  each  of  the  coefficients  an] 
would  need  to  be  doubled  as  well. 


In  addition,  the  new  representation  scheme  includes  a  spline  representation  of  the  design 
waterline.  In  the  old  scheme,  it  was  assumed  that  the  design  waterline  was  the  line  s  =  0; 
however,  this  is  not  practicable  if  one  wishes  to  allow  modifications  in  the  draft  waterline  of 
modified  hulls.  Instead,  as  described  in  Section  8,  the  design  (or  draft)  waterline  is  defined  by 
s  =  w(z'),  where  w(x/ )  is  a  function  approximated  by  a  spline.  The  spline  data  are  included 
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in  the  data  file  immediately  following  the  spline  coefficients  describing  the  hull.  If  a  hull  data 
file  read  by  a  HLLFLO  program  does  not  include  design  waterline  data,  the  HLLFLO  program 
will  assume  that  the  design  waterline  is  s  =  0  and  that  Zw  =  D.  The  format  for  the  design 
waterline  data  is 

Variable 

ZWTR 

KWTR,  NWTR 

TWTR(J),  J=1,...,KWTR+NWTR 
BWTR(J),  J=1,...,NWTR 

where: 

•  ZWTR  is  the  design  draft  of  the  ship, 

•  KWTR  is  the  order  of  the  spline  used  to  define  the  position  of  the  design  waterline, 

•  NWTR  is  the  number  of  B-splines  used  to  define  the  position  of  the  design  waterline, 

•  TWTR  is  the  knot  sequence  used  to  define  the  position  of  the  design  waterline,  and 

•  BWTR  are  the  spline  coefficients  used  to  define  the  position  of  the  design  waterline. 

The  program  SMHULL  has  been  modified  to  allow  the  user  to  spline  the  design  waterline  after 
the  splines  for  the  hull  geometry  have  been  determined.  The  user  prompts  and  sub-commands 
allowed  for  spline  modification  are  exactly  as  described  in  Section  A.  8. 

All  the  HLLFLO  programs  have  been  modified  so  that  they  use  the  new  representation 
scheme.  These  changes  should  be  transparent  to  users  but  for  the  fact  that  the  spline  coefficients 
in  the  hull  data  files  (the  SMHULL  output  file)  will  have  changed  by  a  factor  of  B  or  D  (the 
format  of  the  data  files  remains  the  same).  Hull  data  files  created  using  HLLFLO  version  1.0 
will  need  modification  to  be  used  with  the  updated  programs.  The  program  CNVHLL  has  been 
provided  as  a  means  for  modifying  the  files.  CNVHLL  simply  asks  for  the  name  of  the  old  data 
file  and  modifies  the  spline  coefficients  within  it  so  that  they  obey  the  new  hull  representation 
scheme.  The  new  data  is  written  into  a  file  specified  by  the  user. 

Except  where  noted  in  Chapter  5  of  the  HLLFLO  User’s  Guide[l],  CNVHLL  is  written  in 
ANSI  standard  FORTRAN  77.  It  calls  subroutines  from  the  spline  library  BSPLIN[5],  from  the 
GETWRD  Package(9j,  and  from  the  HLLFLO  libraries  CMDSUB,  HCRSUB,  and  MSGSUB. 
Provided  that  the  PCL  command  HLLFLO  has  been  declared  (see  Appendix  A),  CNVHLL  may 
be  run  on  the  DREA  DEC-20/60  using  the  command 


Format 

(G15.7) 

(215) 

(5G15.7) 

(5G15.7) 


fl HLLFLO  CNVHLL 
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